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Abstract. Measurements of the cross section for exclusive virtual-photoproduction of p mesons from 
hydrogen are reported. The data were collected by the HERMES experiment using 27.5 GeV positrons 
incident on a hydrogen gas target in the HERA storage ring. The invariant mass W of the photon-nucleon 
system ranges from 4.0 to 6.0 GcV, while the negative squared four-momentum Q'^ of the virtual photon 
varies from 0.7 to 5.0 GeV^. The present data together with most of the previous data in the intermediate 
VF-domain are well described by a model that infers the W-dependence of the cross section from the 
dependence on the Bjorken scaling variable x of the unpolarized structure function for deep-inelastic 
scattering. In addition, a model calculation based on Off- Forward Parton Distributions gives a fairly good 
account of the longitudinal component of the (P production cross section for > 2 GeV^. 
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1 Introduction 

This paper presents cross section measurements for exclu- 
sive diffractive production of p'^{770) vector mesons in po- 
sitron scattering on a target. The production of vector 
mesons by real or virtual photons is of considerable inter- 
est, as the corresponding cross section is closely related to 
other observables in lepton scattering. For example, with 
help of the recently introduced Off- Forward Parton Distri- 
butions (OFPDs) , one can relate elastic nucleon form fac- 
tors, deep-inelastic scattering (DIS) structure functions, 
virtual Compton scattering cross sections, and vector me- 
son production cross sections . The OFPDs represent 
a generalization of the parton distributions measured, for 
instance, in inclusive DIS experiments. 

In Refs. the longitudinal part of the virtual- 

photoproduction cross section is calculated in the OFPD 
framework. At large values of the photon-nucleon invari- 
ant mass {W > 10 GeV) the calculated cross section is 



dominated by a two-gluon exchange mechanism, which has 
been treated using the perturbative approach of Ref. [Q . 
These calculations reproduce existing data. However, un- 
certainties arise due to the size of higher-order and higher- 
twist contributions , which are larger at small values 
of (the negative square of the four-momentum of the 
virtual photon). Below 10 GeV the calculated cross section 
is dominated (in leading twist) by a handbag diagram, in 
which the virtual photon is absorbed by a valence quark 
in the nucleon |^,^. Following Ref. Q this mechanism 
is called quark exchange. However, few data are available 
between 4 and 10 GeV. Several data sets exist for W < 
4 GeV, but in this domain the OFPD calculations do not 
apply, as the reaction receives contributions from other 
reaction channels in this domain. 

The relation between the vector meson leptoproduc- 
tion cross section and the structure function -F|'(a;) of the 
proton (with x the Bjorken scaling variable) is apparent 
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in a model calculation by Haakman et al. based on 
Reggeon field theory. In this model the x-dependence of 
F2{x) is related to the W^-dependence of the cross sec- 
tion for vector meson leptoproduction. The model gives a 
good description of j*p — > p^p cross sections for W rang- 
ing from 10 to 180 GeV. Between 4 and 10 GeV the cal- 
culated shape of the W^-dependence changes significantly 
with Q^. The model has not yet been properly assessed in 
this domain due to lack of data. 

The two examples above illustrate the need for addi- 
tional exclusive leptoproduction data in the W^-range 
between 4 and 10 GeV. Such data are reported in the 
present paper. Exclusive cross sections have been mea- 
sured for W values between 4.0 and 6.0 GeV, and 
values between 0.7 and 5.0 GeV^. The kinematical vari- 
ables used to describe leptoproduction are introduced 
in the next section. Experimental details are provided in 
section 3, while the data analysis is discussed in section 4. 
The experimental results are compared to various model 
calculations in section 5, and the paper is summarized in 
the last section. 



2 Kinematics 

In the present leptoproduction measurement both the scat- 
tered lepton and two oppositely charged hadrons are ob- 
served. The leptoproduction events are identified by 
requiring that the reconstructed invariant mass of the two 
hadrons is close to the mass of the p° meson, i.e. 0.77 
GeV. The kinematic variables that characterize such mea- 
surements are defined here in the laboratory frame. The 
variables describing the kinematics of the virtual photon 
include its energy its fractional energy y = v/Eq (with 
Eq the incident lepton energy), and its four- momentum q. 
The latter quantity is related to the previously introduced 
variable through — -q^ > 0. The Bjorken scal- 
ing variable x is related to ly and hy x — Q"^ / {2Mv), 
and the photon-nucleon invariant mass W is given by 
W'^ = M"^ + 2Mv - Q^. The virtual-photon polarization 
parameter is represented by e. 

The p*^ meson is characterized by its four-momentum 
V and ^TTTT) which represents the angle between the two 
pions into which the p° meson decays. 

Combining some of the lepton and meson variables, 
one can introduce t = {q — v)"^ < 0, the square of the four- 
momentum exchange between the virtual photon and the 
target, and AE = (P^ - M'^)/2M, which is a measure of 
the missing energy (where Py ^ P + q — v represents the 
4-momentum of the unobserved final state Y with P = 
(M, 0) that of the target nucleon). 

Instead of —t, the above-threshold momentum transfer 
—t' = —t -\-to is often used, which is approximately equal 
to - the square of the transverse momentum of the p^ 
meson with respect to the direction q. In this expression 
—to represents the minimum value of — i for fixed values 
of ly, and Py- 



3 Experiment 

The data were collected during the 1996 and 1997 running 
periods of the HERMES experiment flltl at DESY using 
a 27.5 GeV longitudinally polarized positron beam with 
a gas target in the HERA storage ring. Part of the 
data set was collected with longitudinally polarized tar- 
gets. Since the polarization degrees of freedom were not 
exploited in the present analysis, the average over both 
target polarization states is taken. 

The HERMES polarized proton target ^ is formed 
by injecting a nuclear-polarized beam of atomic hydrogen 
from an atomic beam source into a tubular open-ended 
storage cell inside the positron ring. The cell provides a 
40 cm long target of pure atomic species with an areal 
density of approximately 7 x 10^^ atoms/cm^. With un- 
polarized targets of molecular hydrogen, areal densities 
of about 10^^ atoms/cm^ were obtained. The storage cell 
was shielded from synchrotron radiation by two sets of 
collimators, one of which is moveable. No particles were 
observed to originate from scattering events in the walls 
of the storage cell. 

During the course of one fill (typically 8 hours long), 
the positron current in the ring decreased from typically 30 
- 40 mA at injection to 10 mA, at which point the ring 
was emptied. The data presented in this paper correspond 
to an integrated luminosity of 108 pb~^. 

The HERMES spectrometer is described in detail else- 
where |ll| . It is a forward spectrometer in which both 
the scattered positron and produced hadrons are detected 
within an angular acceptance ± 170 mrad horizontally, 
and ± (40 - 140) mrad vertically. The scattered-positron 
trigger was formed from a coincidence between a pair of 
scintillator hodoscope planes and a lead-glass calorime- 
ter. The trigger required an energy of more than 3.5 GeV 
deposited in the calorimeter. (For part of the running 
the trigger threshold was reduced to 1.5 GeV, for which 
a correction was applied in the data analysis.) Positron 
identification was accomplished using the calorimeter, the 
preshower counter consisting of the second hodoscope pre- 
ceded by a lead sheet, a transition-radiation detector, and 
a threshold gas Cerenkov counter. This system provided 
positron identification with an average efficiency of 99% 
and a hadron contamination of less than 1%. 



4 Data analysis 

Only those events were selected that contained a scattered 
positron and exactly two hadrons with opposite charge. 
(A more detailed description of the analysis is given in 
Rcfs. [^|l^.) A number of geometric requirements were 
imposed on the particle tracks to ensure that they were 
well contained within the acceptance of the spectrometer. 
It was also required that the tracks originated from along 
the beam line within ± 18 cm of the centre of the target. 
In addition, several constraints were imposed on the kine- 
matic variables. The size of the radiative corrections was 
limited by requiring y < 0.85. Because the VF-acceptance 
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Fig. 1. Missing energy spectra for leptoproduction of me- 
sons off at an incident energy of 27.5 GeV. The upper 
(lower) panel corresponds to the W-range between 4 and 5 
(5 and 6) GeV. The data have been selected requiring 0.6 GeV 
< < 1 GeV and ~t' < 0.4 GeV^. The solid histograms 

represent the results of a Monte Carlo simulation |^ of the 
background, which has been scaled to the data for AE > 3.0 
GeV. The dotted histogram includes a simulation jl4j of exclu- 
sive p'^-production as well. As the Monte Carlo simulations do 
not include inelastic nucleon excitations or internal radiative 
effects for the exclusive channel, the Z\£'-region between 0.5 
and 2.5 GeV cannot be properly described. 



of the HERMES spectrometer for p° production is sharply 
reduced both below 4 GeV and above 6 GeV, cross sec- 
tions for only two M^-bins (4-5 and 5-6 GeV) have been 
extracted from the data. 

The vector mesons were identified by requiring 0.6 
GeV < Mtttt < 1 GeV, with A/tttt the invariant mass of the 
pair of detected hadrons, assuming that they are pions. 
It has been verified that this requirement also removes 
the (f) — > K+K^ background, by confirming that the 4> 
events, which occur at Mkk ~ M^, appear in the Ad-^j^ 
spectrum at Mt^t^ < 0.6 GeV. Here Mkk is the invariant 
mass calculated assuming that the two hadrons are kaons, 
and = 1.019 GeV is the mass of the 4> meson. 

In exclusive p° electroproduction, i.e. eA^ — > e'p'^N, 
the part of the final state that is unobserved at HERMES 
consists of a nucleon recoiling without excitation. Such 
events are selected by requiring that the missing energy 
AE is approximately zero. In this domain the AE spec- 
trum, as displayed in Fig. |l| for the two H^-bins, shows a 
clear peak near AE « 0. The exclusive events were se- 
lected by requiring AE < 0.4 GeV. 



Using the Monte Carlo simulation for exclusive p° pro- 
duction described in section 4.2, it has been evaluated that 
on average 18% of the exclusive events fall outside the im- 
posed AE window. This effect has been accounted for (bin 
by bin) in the acceptance correction. 

The requirement on AE is relatively tight in order to 
suppress inelastic contributions involving nucleon excita- 
tions. As no Monte Carlo simulation is available for such 
double diffractive processes, their remaining contribution 
is estimated from other data in section 4.3. 



4.1 Background subtraction 

There is a background contribution under the exclusive 
peak that is caused by hadrons from DIS fragmentation 
processes. Part of this background is removed by excluding 
the region —t' > 0.4 GeV^, where the background domi- 
nates the p*' yield. The remaining background is estimated 
and subtracted using the methods described below. 

The background estimate is based on the LEPTO Monte 
Carlo program |l^. This simulation includes hadrons re- 
sulting from fragmentation processes in deep- inelastic scat- 
tering (using the Lund fragmentation code), but not from 
diffractive p*^ production. The AE spectrum was gener- 
ated independently for each {Q^ ,W)-hm, and normalized 
to the experimental data in the AE region above 3.0 GeV 
(see Fig. Thus normalized, the background was sub- 
tracted from the data in the AE domain that was used for 
the determination of the cross section. Depending on kine- 
matics the background contribution ranges from (3 ± 1)% 
to (9 ± 4)%, where the uncertainties include the effect of 
varying the zf\_E-ranges used in evaluating the background 
contribution. The lower limit of the Z\i?-range used to nor- 
malize the background was varied from 3 to 5 GeV, and 
the upper limit used in selecting the exclusive events was 
varied from 0.3 to 0.6 GeV. 

The background subtraction procedure was verified by 
also using an alternative method. In this case a AE spec- 
trum was extracted from the data by requiring — t' val- 
ues between 0.7 and 5 GeV^, i.e. well beyond the region 
used to determine the diffractive cross section. A similar 
normalization procedure is used, i.e. the background spec- 
trum is normalized such as to reproduce the same number 
of events for AE > 3 GeV as in the spectrum obtained 
with the standard —t' requirement. The tail of the normal- 
ized background spectrum for AE < 0.4 GeV was taken 
to represent the actual background contribution to the p^ 
mass peak. The cross sections obtained after this back- 
ground subtraction method were found to be consistent 
with those obtained using the Monte-Carlo based method 
within the systematic errors listed below in Table |^. 



4.2 Acceptance correction 

The data were corrected for the finite acceptance and in- 
efficiencies of the HERMES spectrometer for produc- 
tion. The correction was evaluated using a Monte Carlo 
simulation based on the Vector Meson Dominance (VMD) 
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Fig. 2. Invariant mass spectrum of the two oppositely charged 
hadrons assuming that they are pions. Backgrounds have been 
subtracted, the contribution from decays has been removed 
by requiring that Mkk > 1-04 GeV, and the acceptance cor- 
rection has been apphed. The mass peak at 770 MeV is 
described by a Breit-Wigner function (solid line), of which the 
parameters were fitted to the data. 



For each kinematic 
-t') good agreement 



model. Details are given in Ref. | 
variable [i^, Q^, W, A'Ut^, cos^^tt ana 
between the measured distribution and the simulation was 
obtained. 

The acceptance correction is large because the HER- 
MES spectrometer has a relatively small acceptance for 
-f*p p^p events. The acceptance for p° mesons (given 
an observed scattered positron) ranges from only 5.9% at 
low to 18% at high Q^. As the corresponding correction 
factors are large, it is mandatory to study the dependence 
of the acceptance correction on the details of the Monte 
Carlo calculation. For that purpose, the entire analysis 
was redone varying the assumptions for the angular dis- 
tributions of the production and decay of the p'^-meson. 
The original Monte Carlo simulation, which assumed s- 
Channel Hclicity Conservation (SCHC), was replaced by 
one that produces uniform angular distributions for all rel- 
evant angles. In addition, the Q'^-dependence of the VMD 
propagator (as given by Eq. |^ in section 5.1) was varied by 
changing the exponent from 2 to 2.5, and an alternative 
event generator DIPSI ||T^ was used. The sensitivity of 
the Monte Carlo yield to variations of the angular distri- 
butions and changes of the assumed Q^-dependence leads 
to typical systematic uncertainties of 10% and 6%, respec- 
tively. The effect of the event generator itself amounts to 
typically 8%. Including all variations, the uncertainty of 
the cross sections due to the acceptance correction ranges 
from 12% to 17%, depending on the kinematics. 



4.3 Cross section determination 

The Mtt-k mass spectrum is shown in Fig. |^. In this spec- 
trum backgrounds are subtracted, the acceptance correc- 



tion is applied and all kinematic requirements mentioned 
above are imposed. The p° mass peak was fitted using 
a Breit-Wigner (p-wave) function and a constant back- 
ground term representing the non-resonant physics back- 
ground. The Breit-Wigner function was multiplied by the 
factor (Mp/M^7r)"° in which Ug represents the skewing pa- 
rameter. This procedure is identical to the one described 
in Ref. |l^, and was first proposed by Ross and Stodol- 
sky ijl^. The resulting fit parameters (M^o — 0.776 ± 
0.003 GeV, Fpo = 0.147 ± 0.009 GeV, and = 2.2 ± 
0.4) are in good agreement with the values published by 



E665 (T|: Mpo = 0.777 ± 0.002 GeV 



Fpo = 0.146 ± 



0.003 GeV, and n, = 2.7 ± 0.6 (for the same Q^.region). 
Both experiments find a background term that is consis- 
tent with zero, and are in fair agreement with the PDG- 
values |l|: MpO = 0.7700 ± 0.0008 GeV, Fpo = 0.1507 ± 
0.0011 GeV. If the effect of the cu decay is also considered 
in the fit, lu mass and width parameters are found that 
are also consistent with Ref. [|l^ . The contribution of this 
interference averages out to zero in the considered mass 
interval. The absolute contribution of w-production to the 
apparent p° yield is estimated to be less than 1%, and has 
therefore been neglected. It should be noted that the p° 
production cross sections presented in this paper corre- 
spond to the data in the mass region 0.6 GeV < Mt^-^ < 
1.0 GeV with a correction (of 14 ± 2%) for the tails of the 
skewed Breit-Wigner distribution outside this mass 
window. This correction factor has been evaluated rela- 
tive to the mass range between 2M7r and {MpO + Sr'po). 

Examples of the —t' distribution of the data can be 
found in Refs. [|l^,^. The —t' distribution is fit to a falling 
exponential, and the fit is used to correct the measured 
cross section for the excluded region —t' > 0.4 GeV^. The 
slope parameter b from the fit is (6.82 ± 0.15) GeV~^, 
which is consistent with previously published values. 

The production data were normalized using the in- 
clusive deep-inelastic scattering (DIS) yield derived from 
the same data sample. By comparing the DIS yield to the 
LEPTO Monte-Carlo yield based on the world data [|ll 
evaluated at the appropriate and x values, the required 
normalization factors were determined. Virtual photopro- 
duction cross sections, cr^.p^pOp, were obtained from the 
leptoproduction cross section after applying the photon 
flux factor Ft (using the Hand convention as outlined in 
Ref. [EtI, for instance): 



1 da„ 



Ft dvdQ^ 



(1) 



In doublc-diffractive (DD) processes the target nucleon 
is broken up. In a recent study of target dissociation in 
p^ production at center-of-mass energies between 60 and 
180 GeV, the cross section ratio of double-diffractive to 
single-diffractive production was measured to be 0.65 
± 0.17 1^^. Accounting for the shape of the baryonic spec- 
trum (taken from pp diffractive scattering |2^) and the 
resolution of the present experiment, the DD contribution 
to the diffractive production cross section at HER- 
MES is found to be (4 ± 2)%, for which the data were 
corrected. The uncertainty in the DD contribution origi- 
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Table 1. The measured virtual-photoproduction cross sections 
o"^«p^pOp for exclusive production on (in ^b) corrected 
for radiative effects. Both the statistical (first) and systematic 
(second) uncertainties are listed. 



{Ql [GeVi 


{W) = 4.6 GeV 


{W) = 5.4 GeV 


0.83 
1.3 
2.3 
4.0 


2.46 ± 0.13 ± 0.51 
0.92 ± 0.03 ± 0.15 
0.43 ± 0.02 ± 0.08 
0.16 ± 0.01 ± 0.03 


2.04 ± 0.10 ± 0.43 
1.00 ± 0.04 ± 0.16 
0.41 ± 0.02 ± 0.07 
0.10 ± 0.01 ± 0.02 



nates from a rough estimate of the acceptance for double 
diffractive events relative to that for single diffractive 
events, and the uncertainty of the data of Ref. Taking 
into account other estimates of the ratio of double diffrac- 
tive to single diffractive p° production also at lower W- 
values , which are smaller in general, a total systematic 
uncertainty of about 3% (relative to the measured cross 
section) has been assigned to the DD contribution. The 
small size of the DD contribution is related to the good 
energy resolution of the HERMES experiment {aAE ~ 
0.25 GeV - see Fig. |l|). 

The final cross sections are obtained by applying a 
radiative correction to the experimentally determined 
production cross section. The internal radiative correction 
has been evaluated separately for each {Q^,W) bin, and 
typically amounts to 18% [g5[. External radiative effects 
(caused by detector materials) are included in the accep- 
tance Monte Carlo. 



5 Results 

In Table |] the results for the virtual-photoproduction 
cross section are given for each and W bin. The system- 
atic uncertainties are dominated by those from the accep- 
tance correction factors, which amount to 17%, 11%, 14% 
and 17% in the bins centered at values of 0.83, 1.3, 
2.3 and 4.0 GeV^, respectively. Another important con- 
tribution to the systematic uncertainty stems from from 
the uncertainty in the reconstruction and data selection 
efficiency, which totals about 9%. This contribution has 
been estimated by varying the large number of require- 
ments used to select data and reconstruct valid events. 
The quoted 9% also includes uncertainties in tracking effe- 
ciency, kaon contamination, radiative corrections and DD- 
contribution. Remaining contributions to the total sys- 
tematic uncertainty are the absolute normalization (6%) 
and the background subtraction (4%). The combined sys- 
tematic uncertainty on the yO° cross sections ranges from 
16 to 21%, depending on the kinematics. 

In Fig. ^ the virtual-photoproduction cross section is 
plotted versus W, and compared to existing measurements 
[|l7|, ^6|-^ at nearby values of W. As the various data 
sets have been measured at different average values, 
most data sets were rescaled to the average values of 
the HERMES data using the Q^-dependence of the VMD 
model, which is in agreement with the present data as 
will be shown below (see Fig. ^ and Eq. (H) with m = 
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Fig. 3. The virtual-photoproduction cross section for exclu- 
sive production versus W for the indicated values of (Q^). 
The data collected in the present experiment are represented by 
solid circles. The open squares are from Ref. teSK the open tri- 
angles from Ref. [M , the open stars from Ref!~|2j , the crosses 
from Ref. the open circles from Ref. and the open 
diamonds from Ref. pol]. Previous data have been scaled to 
the presently used Q -bins using Eq. (2). The error bars in- 
clude both the statistical and systematic uncertainties (added 
in quadrature). 



2). It is noted that the various data sets have also been 
obtained at different average e- values, for which no correc- 
tion has been applied. The HERMES data are seen to fill 
the gap which previously existed between the steeply de- 
clining data at low W, and the data collected at higher W 
values, which have a much flatter VF-dependence. With 
a few exceptions, there is a fair consistency among the 
world's data on exclusive p'^ production. The first excep- 
tion concerns the highest W point (at 3.5 GeV) of the 
Cornell data H at = 0.83 GeV^, which deviates by 
about 2(7-3(7 from a smooth interpolation of all other data 
in that M^-region. Secondly, there is a difference of about 
a factor of two between the NMC Q and E665 data 
at ~ 15 GeV in the two highest Q^-bins. The latter 
discrepancy has been reported before |1^, |T| . Part of the 
discrepancy might be related to a (model-dependent) sub- 
traction of the non-resonant contribution to the p°-peak 
that was applied by NMC ||o| in contrast to E665 Q 
and the present experiment. 

It should be noted that it is the total cross section 
for the process ep — > eh+h^p with Mt^tt ~ MpO that is 
displayed in fig. || for all data sets. This total cross sec- 
tion is believed to receive contributions from exclusive p*^ 
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Table 2. The fitted pliotoproduction cross sections for pro- 
duction on ^H. In tlie top row tlie results are sfiown witfi tlie 
variables ao and of Eq. (2) as free parameters, the middle 
row lists the results obtained with ctq and m left free, and the 
bottom row shows the fitted values of m if a parameterization 
of the world data is used to constrain ao The listed errors in- 
clude both the statistical and systematic uncertainties, added 
in quadrature. In all cases the x^/d-o.f. values of the fits are 
close to or less than unity. 



Variables 


{W) = 4.6 GeV 


(W) = 5.4 GeV 


(To [fih] 


10.7 ± 1.8 


13.1 ± 1.7 


e 


-0.023 ± 0.049 


-0.13 ± 0.03 


o-Q [/xb] 


11.2 ± 2.9 


19.7 ± 4.9 


m 


2.42 ± 0.18 


2.83 ± 0.18 


(To [^b] 


11.5 (fixed) 


11.0 (fixed) 


m 


2.44 ± 0.05 


2.46 ± 0.06 



production through Reggeon exchange, as well as from 
reaction channels such as those involvin g n ucleon reso- 
nances at low W, as reporte(i in Refs. Pq-psf. In virtual- 
photoproduction it is usually assumed that the latter con- 
tributions are negligible beyond W ~ 4 GeV. The data 
of Fig. ^ indicate that these additional processes have a 
decreasing contribution to p° production up to about 4-5 
GeV, where the steep decrease of the cross section changes 
into an almost flat VF-dependence. Henceforth, we restrict 
the comparison with existing calculations to data collected 
atW > A GeV. 



5.1 Q^-dependence 

The HERMES data are displayed in Fig. || as a function of 
for the two W bins. For the purpose of extrapolating 
to — 0, the Q^-dependence of the data has been pa- 
rameterized using the following functional form, inspired 
by the VMD model il|,|| 



a(Q') = do 



Af2 



{l + eR{Q^)). (2) 



The VMD-model predicts that m = 2 and that the ra- 
tio of longitudinal to transverse p*' photoproduction cross 
sections is given by i? = ^^Q^/M^. A good description of 
the data is obtained by fixing m to its VMD value, and 
treating ag and as free parameters. Here ag represents 
the p° production cross section for real (Q^ = 0) photons. 
The results of the fit (solid curves in Fig. ^) are listed in 
table ^. 

As the ratio R must be positive, the negative fit val- 
ues of indicate a deficiency of the VMD-model. Hence, 
a different approach was also used, i.e. the exponent m 
of the propagator in Eq. (0) was used as a free parame- 
ter, while the value of R{Q^) was fixed using information 
extracted from measurements of the decay angular dis- 
tribution fill: 



10" 



o 

a. 



10-1 



10° 



10 



1 



<\\ > = 4.6 GeV 



<W> = 5.4 GeV 




R{W,Q^) = c^iW) ■ {Q^ /MlfK 



(3) 



Q2 [GeV^] 

Fig. 4. The virtual-photoproduction cross section for p° 
production on versus for 4 GeV <W < h GeV (upper 
panel) and 5 GeV <W < Q GeV (lower panel). The error bars 
include both the statistical and systematic uncertainties, added 
in quadrature. The solid, dashed and dotted curves represent 
fits using the values of ctq and ao and m, or only m as free 
parameters, respectively. 



The parameterization of R expressed by Eq. (^) has been 
obtained from a fit of the world data on i?, yielding Ci 
= 0.61 ± 0.04, and cq{W) = 0.33 ± 0.03 (0.48 ± 0.03) 
for 4 GeV <W <7 GeV (W^ > 7 GeV). The resulting 
fits (dashed curves in Fig. H) also give a good account of 
the Q^-dependence of the p production cross section. The 
fitted parameters are listed in the middle row of table |^. 
The values of (Tq are somewhat different from the results 
of the previous method, but are consistent as the errors 
have increased considerably. 

It is concluded that an increase of the exponent m 
can avoid the negative values of found in the previous 
fits. This result quantitatively confirms similar conclusions 
obtained by the E665 collaboration (m = 2.51 ± 0.07 at 

« 15 GeV (TH). 

The fitted values of the extrapolated photoproduction 
cross section (Tq are in agreement with existing real-photon 
data ]3^ , but carry larger error bars particularly if the 
exponent and the cross section are treated as free param- 
eters. Since the VMD cross section scales with Q"^ while 
a perturbative QCD description predicts that the lon- 
gitudinal part of the cross section scales with , it is 
of particular interest to study the value of the exponent 
of the propagator in the transitional domain probed by 
the present data. By fixing the value of the photoproduc- 
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W [GeV] 

Fig. 5. The virtual-photoproduction cross section for pro- 
duction versus W at average values of 0.83, 1.3, 2.3 and 4.0 
GeV^ (from top to bottom). The coding of the data symbols 
is the same as in Fig. 2. The solid (and dashed) lines represent 
the results of calculations by L. Haakman et al. For the 
solid (dashed) hues the HERMES and E665 (NMC) data were 
used to normalize the curves. Note that both the data and the 
calculations have been multiplied by the factors indicated on 
the left-hand side of the figure for plotting purposes only. 



tion cross section ctq using a parameterization of the exist- 
ing photoproduction data |Q , a one-parameter fit of the 
Q^-dependence of the cross section has been carried out, 
yielding fairly precise values for the exponent m as can 
be seen from the bottom row of table The fitted value 
is remarkably close to 2^ which is in agreement with re- 
sults obtained by E665 fT^] (2.51 ± 0.07) at an average W 
value of 17 GeV, but somewhat above the result obtained 
by HI H (2.24 ± 0.09) at = 75 GeV. 

5.2 ly-dependence 

In Fig. ^ the virtual-photoproduction cross sections are 
shown as a function of W for the four different bins. As 
was argued above, only the data for > 4 GeV are being 
considered. The data are compared to the calculations of 
Haakman et al. , which relate the M^-dependence of the 
vector-meson production cross section to the proton struc- 
ture function _F|'(a;,Q^). In order to understand how the 
a;-dependence of F|'(x, Q^) influences the M^-dependence 
of the p° production cross section, it is useful to realize 
that X and W are related by W"^ ^ Q'^^ + M"^, thus 
showing that the two variables are essentially inversely 



proportional at large values of and small values of x. 
The relation between the p° production cross section and 
the structure function is apparent in the expression for the 
f-dependence of the cross section given in Ref. [|lo| : 

^ = fp{Q^)F^{x,Qyexp[A{W)t]. (4) 

In this expression the scaling function fp{Q^) is not deter- 
mined by the model, and has been obtained from a fit to 
the experimental data displayed in Fig. |[ F or Fj'(a;, Q^) 
the parameterization of Capella et al. ||3q ] in terms of 
Pomeron and Reggeon exchanges has been taken. The 
slope parameter A is given by 

A{W) = 2[Rl^ + 2a'p \n(W/Wo)] (5) 

with the radius parameter R^^ = Rp + Ry{Q^) expressing 
the combined size of the scattering objects, and Wq given 
by Af^o +Q^. Expressions for a'p, the slope of the effective 
Pomeron trajectory, and Ry{Q^) are given in Ref. [p^ . 
Thus evaluated, the value of A turns out to be consistent 
with the value of the slope parameter derived from a fit 
of the — f distributions | p^ . 

In order to compare the VF-dependence of the model 
of Ref. to that of the data, the calculations have been 
normalized to the HERMES and E665 data, for each 
bin separately. The results are displayed in Fig. || for 
bins centered at 0.83, 1.3, 2.3 and 4.0 GeV^, respectively. 
A good description of these data is obtained, thus confirm- 
ing the findings of Ref. ||l^ but now at lower values of W. 
It should be noted, however, that the data collected at W 
< 4 GeV - as displayed in Fig. || - show a steep rise with 
decreasing W that is not reproduced by the calculations 
of Ref. |lC|]. As was mentioned above, this is presumably 
caused by additional reaction processes that are not con- 
tained in the theoretical framework of Ref. jlO| . Without 
explicit calculations that include these effects, it cannot 
be excluded that there is a finite contribution due to such 
processes above W — A GeV, especially at the lowest 
values. 

The NMC data (also shown in Fig. |) can be com- 
pared to the calculations of Ref. [0 as well. This leads to 
differently normalized curves (dashed lines in Fig. If 
the NMC normalization is adopted, the curves fall below 
the E665 and HERMES data. Within the constraints of 
the present model calculation, the E665 normalization is 
therefore preferred. However, on the basis of the HERMES 
data alone no such distinction can be made. 



5.3 The longitudinal cross section 

The results of recent OFPD calculations, such as those de- 
scribed in Refs. |^-|^, only concern the longitudinal com- 
ponent (Ji of the p'^ virtual-photoproduction cross section, 
because the factorization theorem |^ applies only to the 
longitudinal case. Hence, values oIul were extracted from 
the data. 
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Fig. 6. The longitudinal component of the virtual- 
photoproduction cross section for production versus W at 
average values of 2.3 (left) and 4.0 GeV^ (right). The solid 
hues represent the results of the calculations of Refs. ^ . 
The dashed (dotted) curves represent the quark (two-gluon) 
exchange contributions within these calculations. 



The longitudinal cross section (Tl is related to the total 
p° production cross section aj*p_,pOp: 



R 



1 + eR 



'7*p — >p^'p • 



(6) 



Assuming SCHC the ratio R — ctlI^^t can be derived 
from the longitudinal fraction Tqq of mesons, which can 
be extracted from their polar decay angular distribution 
W(cos 0). In Ref. |l4j such an analysis is presented, which 
has provided us with the parameterization of R shown in 
Eq. (|^). This has been used to evaluate aj^ for each of 
the cross sections listed in Table |l|. (Note that e can be 
determined from the kinematics of each event.) The results 
are listed in Table || and, for > 2 GeV^, displayed in 
Fig. ^ It may be noted that the values of gl increase by 
about 30% if the quoted parameterization of R for W > 
7 GeV is used. 

The data are compared to the model calculations de- 
scribed in Ref. Q , but with the improvements described in 
Ref. [|j : higher twist effects are included in a phenomeno- 
logical fashion, the Q^-dependence of the strong coupling 
constant as was accounted for, and the sea quarks are 
included using the MRST98 parton distributions |3^] . As 
can be seen from Fig. 17 of Ref. ^ the higher-twist con- 
tributions, which were included through the dependence 
on the intrinsic transverse momentum k±^ , reduce the lon- 
gitudinal cross section by typically a factor of 5 at « 
2 GeV^. As the (phenomenological) higher twist correc- 
tions at lower values are even larger, the comparison 
between the OFPD calculations of Ref. |^,|| and the data 
has not been extended to the data below k, 2 GeV^. 
Moreover, the OFPD calculations cannot be used for data 
collected at < 4 GeV, as no other reaction channels 
were included. 



Table 3. Longitudinal cross sections gl for virtual photo- 
production on in pb. The listed uncertainties include both 
the total error on the measured virtual-photoproduction 
cross sections and the error on the parameterization of R for 
W < 7 GeV, which was used in the evaluation of Eq. (^). The 
average values of e for each bin are listed as well. 



< 0" > 


<W> ^4.6 GeV 


<W > ^5.4 GeV 


[GeV^] 


< e > 




< e > 




0.83 


0.87 


0.77 ± 0.16 


0.74 


0.67 ± 0.14 


1.3 


0.87 


0.35 ± 0.07 


0.73 


0.41 ± 0.08 


2.3 


0.86 


0.21 ± 0.04 


0.71 


0.21 ± 0.04 


4.0 


0.83 


0.09 ± 0.02 


0.69 


0.06 ± 0.01 



The OFPD calculations are compared in Fig. ^ to the 
longitudinal cross sections measured by E665 [|7| and the 
present experiment. The calculated contributions due to 
quark exchange (dashed curves) and two-gluon exchange 
(dotted curves) are also shown separately. In the present 
kinematic domain the quark contributions dominate, and 
only at = 4.0 GeV^ and W > 10 GeV the gluon con- 
tribution starts to contribute significantly. A fairly good 
agreement between the calculations and the data is ob- 
tained, given the existing theoretical uncertainties related 
to the size of the higher-twist contributions, and the rel- 
atively low values involved. The calculated rise of the 
cross section at « 5 GeV, which is associated with the 
contribution due to the exchange of valence quarks, is not 
inconsistent with the new HERMES data. 



6 Summary 

In summary, cross sections have been presented for ex- 
clusive diffractive virtual photoproduction in the W- 
domain between 4 and 6 GeV. The Q^-dependence of the 
cross section is well described by the propagator of the 
Vector Meson Dominance model, although an increase of 
the exponent from its original value of 2 to about 2.5 
is needed when consistency with existing data for R — 
ol/<^t is required. By extrapolating the 7*p — > p^p cross 
section to = 0, photoproduction cross sections are 
found which are in fair agreement with previously pub- 
lished values of about 11 ph. The VF-dependence of the 
present and (most) existing data in the VF-domain be- 
tween 4 and 25 GeV is well described by a model link- 
ing the x-dependence of the proton structure function 
F^iXjQ^) to the vector meson production cross section. 
The longitudinal component of the cross section has been 
compared to calculations based on the Off-Forward Par- 
ton Distribution framework. A fairly good agreement with 
the data is found. 
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